Delayed double ionization as a signature of Hamiltonian chaos 
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We analyze the dynamical processes behind delayed double ionization of atoms subjected to strong 
laser pulses. Using reduced models, we show that these processes are a signature of Hamiltonian 
chaos which results from the competition between the laser field and the Coulomb attraction to the 
nucleus. In particular, we exhibit the paramount role of the unstable manifold of selected periodic 
orbits which lead to a delay in these double ionizations. Among delayed double ionizations, we 
consider the case of "Recollision Excitation with Subsequent Ionization" (RESI) and, as a hallmark 
of this mechanism, we predict oscillations in the ratio of RESI to double ionization yields versus laser 
intensity. We discuss the significance of the dimensionality of the reduced models for the analysis 
of the dynamical processes behind delayed double ionization. 
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I. INTRODUCTION 

The dynamical processes by which atoms lose their 
electrons when subjected to strong laser pulses usually 
fall into two main categories [1]: sequential double ion- 
ization (SDI) in which electrons are removed one after 
the other, and a less trivial class of processes coined 
nonsequential double ionization (NSDI). In NSDI pro- 
cesses, the Coulomb interaction between two electrons 
plays an essential role. A particularly relevant example 
is afforded by the recollision scenario in NSDI Q in 
which a preionized electron returns to the ion core to 
excite a bound electron. This excitation may lead to a 
double ionization if the energy transfer from the preion- 
ized electron to the bound electron is sufficient to over- 
come the Coulomb attraction to the nucleus, while the 
preionized electron keeps enough kinetic energy to remain 
ionized. This so-called "direct impact ionization" [4] is 
often invoked to explain the additional amount of dou- 
ble ionization at intermediate values of laser intensities. 
However the dynamics exhibits a much richer variety of 
double ionization processes. It has been observed that, 
sometimes, the recollision excites the parent ion which is 
ionized by the field after some delay, ranging from a quar- 
ter up to several laser cycles. These processes mimic a 
sequential double ionization process even though they fall 
into the NSDI category since the electron-electron inter- 
action plays a role. This alternative and less straightfor- 
ward road to NSDI is called Recollision Excitation with 
Subsequent Ionization (or RESI for short 0-0 )• These 
processes are by no means rare and unimportant. For 
instance, it was shown that RESI plays a central role in 
High Harmonic Generation [8| which is a way to gen- 
erate short X-ray pulses. RESI processes belong to a 
category of double ionization involving a delay in the 
second ionization, and this category is referred to as de- 
layed double ionization. The precise mechanism respon- 
sible for delayed double ionization, is still a subject of 
study from the theoretical and experimental points of 
view [il . Experimentally, the signature of RESI can be 



found in momentum distributions: The two lateral peaks 
in the bimodal distribution are attributed to direct im- 
pact ionization whereas RESI contributes to the central 
part around zero ion momentum 

[J,i[iQHil. It also 
turns out that direct impact ionization is expected to be 
triggered when the electric field is close to zero whereas 
RESI is expected when the electric field is maximum. 

There is still some debate on the mechanisms of de- 
layed double ionization @,|l3|. The conventional scenario 
for RESI is based on a tunneling argument in the adia- 
batic approximation. This scenario is depicted in Fig. [TJ 
After a recollision at a zero of the field, the remaining 
ion is in an excited state but not yet ionized. When the 
laser reaches an extremum, it creates a potential bar- 
rier through which the electron can tunnel. Therefore 
it raises a natural question if RESI is a purely quantum 
phenomenon. To counter this argument, RESI has been 
observed with classical calculations [14-[l6| , for which all 
ionizations are over- 1 he-barrier, casting some doubts on 
the validity of the tunneling argument in the adiabatic 
approximation. 

In this article we revisit this conventional picture of 
the RESI mechanism by analyzing the phase space mech- 
anisms behind delayed double ionizations, of which RESI 
is a prime example. We show that after an electron is pro- 
moted to a specific region of phase space, it is trapped 
for some time before ionizing. This specific region is a 
chaotic one resulting from the competition (and reso- 
nances) of the Coulomb attraction with the laser field. 
At a given laser intensity, this thin chaotic region is or- 
ganized by a single resonance through the hyperbolic pe- 
riodic orbit it creates. The chaotic nature of RESI tra- 
jectories was already identified in Ref. [8] through the 
"deflection function". The way an electron is ionized in 
the chaotic region is by following the unstable manifold 
of some selected periodic orbits. As a consequence, the 
properties of these periodic orbits (as a function of the 
laser intensity) drive the delayed double ionization yields. 
As a signature of these periodic orbits, we find oscillations 
in the RESI to double ionization yields as a function of 
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FIG. 1. (color online) Schematic representation of the con- 
ventional mechanism for RESI using a tunneling argument in 
the adiabatic approximation. Left panel: After a recollision 
with a preionized electron at a zero of the field, the inner elec- 
tron is excited, but not ionized. Right panel: Later on, at an 
extremum of the field, the potential barrier is lowered by the 
laser field, and the excited electron tunnels. On each panel, 
the curves represent the effective potential (soft Coulomb po- 
tential plus the laser) , the dark lines label the laser excitation 
and the horizontal lines represent the energy levels for the 
electron. 

the laser intensity in a range of laser frequency. Some of 
these results were announced in a recent Letter p^Tf . 

In this paper, we restrict our analysis to linearly polar- 
ized fields. Interesting experimental and early theoretical 
results suggest the importance of delayed double ioniza- 
tion with elliptic or circular polarization as well [18]. 

In Sec. ini we describe the classical models we use for 
the analysis of delayed double ionization. We exhibit 



where = {xi^yi^Zi) is the position of the i^^ electron 
{i = 1,2) and = {Px,i^Py,i^Pz,i) its canonically conju- 
gate momentum. Here x and p belong to for (i = 1, 2 
or 3. The nucleus is assumed to be fixed at the origin be- 
cause of its large mass compared to that of the electrons 
and the relative short duration of the laser pulse. The 
linearly polarized laser field (along the direction e^^) is 
characterized by its amplitude Eq and has a wavelength 
of 780 {uj = 0.0584) or 460 nm {u = 0.1 a.u.) which is 
specified for each computation. Its envelope / {t) is com- 
posed of a two laser cycle linear ramp up and a six laser 
cycle constant plateau. The constants a and b are the 
electron-nucleus and electron-electron softening parame- 
ters respectively and are chosen so as to be compatible 
with the ground state energy Eg^ defined as the sum of 
the first two ionization potentials [19, 23-25]. The same 
softening parameters are used independently of the di- 



some relevant reduced models, and show that they are 
pivotal for the identification of the delayed double ion- 
ization mechanism. In Sec. IIII) we present some peculiar 
features of RESI trajectories in phase space. Using the 
reduced models of Sec. [Ill we provide a detailed analysis 
of the mechanisms in phase space. We show the impor- 
tance of the dimensionality of the reduced models. 



II. CLASSICAL MODELS 

Remarkably, classical mechanics is able to identify 
both sequential and nonsequential pathways leading to 
double ionization and reproduces the experimental and 
computational observations closely |19l-[21|. This suc- 
cess was ascribed to the paramount role of the electron- 
electron interaction which is present both in quantum 
and classical models. One of the advantages of classical 
mechanics is to analyze these systems in phase space in 
a very detailed way. In Refs. [21|, i22|] it was shown how 
invariant structures like periodic orbits rule double ion- 
ization processes in strong and short laser pulses. The 
methodology here is to follow the same modus operandi 
in order to provide insights into delayed double ioniza- 
tion. In this section we detail the various models which 
will be used later in the article. 

The parent model is a two active electron atom with 
soft- Coulomb potentials subjected to an intense and 
short Hnearly polarized laser pulse in the dipole approx- 
imation |19l-[21|. The corresponding Hamiltonian reads: 



(1) 

I 

mensionality d for a given atom. 

The model ([T]) has two and a half degrees of freedom. 
In order to describe the dynamical processes in phase 
space, it is desirable to consider reduced models. These 
reduced models are described in what follows. We con- 
sider several kinds of reduced models: dimensionally re- 
duced ones, and single electron models, or a combination 
of both. All of these models are used in a complementary 
way so as to get insights into the dynamical processes of 
delayed double ionizations. 



A. Dimensionally reduced two electron models 

From the dynamics associated with Hamiltonian ([1]), 
we notice that any plane containing the nucleus and the 
polarization axis is dynamically invariant. More specifi- 
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cally, starting with the initial condition zi = Z2 = and 
Pz,i = Pz,2 = 0, the two electrons remain in the (x^y) 
plane. It is then possible to consider a dimensionally 
reduced model where the z-direction (and its canonical 
momentum) is left out. Such models are generally re- 
ferred to as 2 + 2 dimensional models, making reference 
to the dimension of the configuration space in which each 
electron moves. 

The 2 + 2 dimensional model can further be reduced 



to a 1 + 1 dimensional model by constraining each elec- 
tron to move along a line which has a constant angle 
with the axis of polarization. We notice that, contrary 
to the reduction from 3 + 3 to 2 + 2 dimensions which 
can be carried out independently of the ellipticity of the 
field, the one from 2 + 2 to 1 + 1 can be considered only 
for linear polarization. The resulting constrained dynam- 
ics is obtained using, e^., Dirac's theory of constrained 
Hamiltonian systems [26| and the resulting Hamiltonian 
reads 



2 2 

%2 {xi,X2,Pi,P2,t) = y + y 



V4 + V4 + ^(Xi -X2)2+4xi 

+ (xi cosai + X2 cos 0^2) Eq/ (t) sincjt. 
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where is the angle with the axis of polarization (see 
Fig. [2]). Two special cases have been studied in the lit- 
erature. They correspond to ai = 0^2 = 0, where the 
electrons are aligned with the axis of polarization [l9|, [23I - 
[25], and ai = —0L2 = 7r/6 [l^, [28|. In the first case, the 
reduced Hamiltonian reads: 



Hii (x\,X2,P\,P2,'t) 



p\ pI 2 _ 2 

2 2 + a2 ^xl + a2 

1 



^J{x\ - ^2)2 + 
+ {x\ + X2) Eof{t) sincjt. 



(2) 



We note that the reduced dynamics corresponds to a 
symmetry of original 2 + 2 model. In the second case, 
the reduced Hamiltonian reads: 



l-Llie {xi,X2,Pl,P2,t) = 



pI 



pI 



+ {x\ + X2) Eof{t) sincjt, (3) 

where Eq = £^0 a/3/2 = £^ocos±7r/6 is the scaled peak 
field amplitude. Here the lines along which each elec- 
tron is moving are not conserved by the fiow of Hamil- 
tonian ([1]). The main reason why Hamiltonian (|3]) has 
been considered in the literature is to allow the escape 
of electrons simultaneously with symmetrical momenta. 
In the numerical calculations, we use the two models in- 
dependently and show that the dynamical mechanisms 
that regulate delayed double ionization are qualitatively 
the same for both 1 + 1 dimensional models. 

Although the 1 + 1 dimensional reduction may seem 
as a crude approximation, it has shown to be surpris- 
ingly accurate in reproducing quantum and experimen- 
tal results in double ionization qualitatively, such as the 
characteristic knee shape in double ionization yield as a 




FIG. 2. (color online) Dimensional reduction for two active 
electron models. From the 2+2 dimensional model, we con- 
strain the electrons to move along lines that form a constant 
angle with the polarization axis Qx- 



function of the laser intensity. This is the simplest model 
looking at the correlated dynamics of two electrons. In 
this case, the role of the softening parameters is to allow 
the charged particles to pass through each other along the 
polarization axis, as they would have done in a transverse 
direction with a higher dimensional model. 



B. Reduced one electron models 

Although the two electrons are treated equally in the 
previous models, a careful analysis of the dynamics re- 
veals that, given some initial conditions, their dynamics 
are qualitatively very different: typically one electron is 
quickly ionized (the outer electron) while the other one 
remains close to the core (the inner electron) [U, l22j . 
This observation stands independently of the chosen di- 
mension, and it is at the core of the three step sce- 
nario [2, 3]. As a result, for each electron an effective 
reduced model can be crafted. When this outer electron 
is far away from the core the electron-electron interaction 
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can be neglected. Therefore the effective Hamiltonian 
for the inner electron is built from Hamiltonians ([1]), ([2]) 
and dS]): 



T^in (x,p,t) = 




-^xEof (t) sin(a;t+(/)o), 



where Eq = Eq cos a is the scaled laser amplitude in 
1 + 1 dimensional models ([2]) and ([3|), and Eq = 
for higher dimensional models. We notice that, up to 
a scaling of the field amplitude, the same models for the 
inner electron stand independently of the chosen angle 
a. To simplify the comparison between the models, in 
all numerical results we drop the laser amplitude scaling 
parameter (i.e., we choose Eq = Eq). This model for the 
inner electron is only valid after the outer electron has 
left the nucleus (i.e., after the preionization at t = to 
such that mod(to, 27r) = and outside of the recollision 
events). For convenience, we shift time by to so that the 
inner electron model is valid for t > 0. Since most of the 
nonsequential double ionization processes happen in the 
plateau [21], the model for the inner electron is further 
simplified as 

1 1 1 1 2 2 

"Hin (x, p, t) = — , + xEo sm{ujt + (po), 

Vl|xf + a2 

(4) 

where <po denotes the laser phase at which the inner elec- 
tron model is started (e.g., after the preionization to or 
after the end of the final recollision). In what follows, 
we use these reduced models to highlight the relevant 
mechanisms at play for the delayed double ionization. 

III. MECHANISMS FOR DELAYED DOUBLE 
IONIZATION 

A. Delayed double ionization trajectories 

First we begin by analyzing some sample trajectories 
of the two electron Hamiltonian ([1]). For the display, we 
only consider the case d = 2. Other dimensional models 
lead to a similar typology of trajectories and the main 
differences between the various dimensional models are 
analyzed in Sec. lIIIDl Among all the two-electron trajec- 
tories, we consider the ones which are doubly ionized at 
the end of the laser pulse, and which present a significant 
delay between the last time the preionized electron influ- 
ences the core dynamics and the subsequent ionization. 
These trajectories will be referred to as delayed double 
ionizations. We use an ener gy c riterion to identify if an 
electron has ionized or not |2TI , I22I , 0| . The energy £ 
of an electron is defined as the sum of its kinetic energy 
plus the (soft) Coulomb interaction with the nucleus: 




An electron is considered ionized whenever its energy is 
positive, i.e., £ > 0. As for the electron-electron interac- 
tion, there are many ways a preionized electron influences 
the dynamics of the inner electron. The most famous one 
in strong field physics is the recollision. A recollision is 
said to have occurred whenever the distance between the 
two electrons is smaller (or equivalently, whenever the 
Coulomb interaction between the two electrons is larger) 
than some threshold. Of course the definition of the 
threshold is somewhat arbitrary. However on inspection 
of a large sample of trajectories, it leads to the definition 
of a recollision where there is a significant exchange of 
energy between the two electrons. 

We identify three main types of delayed double ioniza- 
tions (to be specified later on) and give an illustration of 
them in Fig. [3l and compare them to the standard direct 
impact mechanism (top left panel). The common fea- 
ture of all the delayed double ionizations shown in Fig. [3] 
is that the inner electron's role is very far from being a 
passive one, waiting for a recollision to ionize as in the 
standard picture of NSDI. In what follows we show that 
the inner electron is promoted to a specific chaotic re- 
gion of phase space where it is trapped for some time 
(which might last up to several laser cycles). The way 
the inner electron is promoted to this region varies ac- 
cording to the specific trajectory and to the values of the 
laser parameters. It is a combination of several mecha- 
nisms: recollision, soft interaction and chaotic diffusion. 
The differences between these three types may seem very 
subtle at first glance. However, by looking closely at the 
trajectory of the inner electron right before being hit by 
the outer electron, we notice that the transfer of energy is 
very low for the soft interaction (lower right panel) while 
it is significant for the recollision (upper right panel). 
The chaotic diffusion (lower left panel) corresponds to 
the extreme case where there is no recollision during the 
plateau (note that during the brief excursion of the firstly 
ionized electron around 1 laser cycles, the electron is not 
ionized). It means that the inner electron is promoted to 
the chaotic layer by the correlated dynamics of the two 
electrons during the ramp up, and no recollision is in- 
volved in the subsequent ionization (leaving out the "R" 
or RESI). Note that chaotic diffusion is not a marginal 
phenomenon among delayed ionization. All three types 
of delayed double ionization share the common featured 
that the inner electron is trapped in a chaotic region of 
phase space for some time before it ionizes. The differ- 
ence between them is the way it reaches this region, i.e., 
through a (soft or strong) recollision or the correlated 
dynamics of the two electrons during the ramp up of the 
field. Among the three types, RESI trajectories, where 
a recollision is involved before ionization, are the most 
famous one. In this case, the significant time delay is 
compared to direct impact double ionization, where the 
subsequent ionization immediately follows the recollision 
(see the upper panels in Fig. [3]). This delay might last 
between a quarter up to several laser cycles. As expected, 
a typical delayed double ionization is the result of a co- 
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FIG. 3. (color online) Typical nonsequential double ioniza- 
tion trajectories for Hamiltonian ([Tj) with d — 2. The pa- 
rameters of the atom are a — h — 1^ and the initial conditions 
have a microcanonical distribution on the energy surface 8g = 
-2.3 a.u. [21^. The laser parameters are / = 3 x lO^^W cm"^ 
and uj — 0.0584 a.u. The top left panel corresponds to the 
conventional direct impact scenario while the other panels 
display the alternative delayed ionization routes (as indicated 
on the figure). Continuous (dashed) lines correspond to the 
dynamics in the polarization (transverse) direction. 



operation between the three mechanisms, so even though 
we have tried to classify these trajectories into three cat- 
egories, the boundaries between them are not sharp at 
all. 



B. Poincare sections 

In order to get deeper insights into the chaotic re- 
gion in which the inner electron is promoted, we consider 
Poincare sections of RESI trajectories of Hamiltonian ([1]) 
with d=l. We have chosen RESI because it constitutes 
a significant subset of delayed double ionization, and it 
is easier to detect numerically. In Fig. |4l we display the 
positions of the electron which ionizes last at the succes- 
sive maxima of the laser field after the last recollision. 
Looking at Fig. IH we note that RESI trajectories are 
restricted to specific parts of phase space. In particu- 
lar, we note a central region, close to the nucleus where 
there is no record of RESI. This region is delimited by a 
thin layer where most of the inner electrons are promoted 
right after a recollision. In addition, we observe swirling 
patterns that extend far away from the nucleus. 

After a recollision, the inner electron is left with the nu- 
cleus and field interactions. We use Hamiltonian (HI for a 
one-dimensional one-electron model. In the left panels of 
Fig. [5l we represent Poincare sections (stroboscopic plots 
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FIG. 4. (color online) Poincare section (stroboscopic plot at 
the maxima of the field) of RESI trajectories after the last 
recollision of Hamiltonian ([2|) with d— 1 . The parameters for 
the atom and the laser are the same as in Fig. O For each 
trajectory, the first point on the section (so right after the 
last recollision) is plotted in light gray (blue online) while the 
following ones are in dark gray (blue). 

at the maxima of the laser, i.e., (j) = ujt -\- mod 27r = 
7r/2) of all the delayed ionizations (i.e., ionizations which 
take more than one laser cycle). For clarity we left out 
the two first points on the Poincare section for each tra- 
jectory (the missing points are located in the same part of 
phase space as the remaining ones, but they clutter the 
fine structures depicted on the figure). Note the close 
similarity between the Poincare section of trajectories of 
the two-electron model (|2]) in Fig. 2] and the ones of the 
one-electron model (j4]) in Fig. [5] (upper left panel). In 
particular the swirling patterns in both figures are nearly 
identical. 

As shown in Figs. S] and [5l RESIs are governed by the 
inner electron model (j4]) after the last recollision. More 
specifically, for one-dimensional models they correspond 
to a thin chaotic layer resulting from a resonance between 
the Coulomb attraction and the laser excitation. Since 
the one-dimensional one-electron model (|4]) captures the 
essential features of delayed ionization, this model is 
studied in more details in the next section. The moti- 
vation for studying one-dimensional models is twofold: 
First, the laser field drives the ionization dynamics along 
the polarization axis. Second, we find that this simplified 
dynamics forms the skeleton of higher-dimensional dy- 
namics which, however, differs from the one-dimensional 
dynamics in significant ways, to be specified in Sec. lIII Dl 

C. One-dimensional one-electron dynamics 

In this section, we analyze the inner electron model 
given by Hamiltonian ^ with only one spatial dimen- 
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FIG. 5. (color online) Left panels: Poincare section (stroboscopic plots at the maxima of the laser field, i.e., = 7r/2) of RESI 
trajectories of the Hamiltonian ^ with d = 1 detected from a large set of random initial conditions. Right panels: Poincare 
section of the stable [light gray (orange online) dots] and unstable [dark gray (blue) dots] manifolds of the periodic orbit 
of Hamiltonian Q with d = 1. The position of the periodic orbit is indicated with a diamond (see the arrow) in the upper 
right panel. The light gray (pink) area in the upper panels represents the part of phase space from which the inner electron 
does not ionize. The lower panels focus are insets of the upper panels. The parameters for the atom and the laser are the same 
as in Fig. [S] 



sion. It has one and a half degrees of freedom with no con- 
served quantity in the presence of a laser field {Eq ^ 0). 
Chaos is expected to develop, but only locally in phase 
space, depending on the parameter values of the laser 
field. To get a rough idea of the dynamics of Hamilto- 
nian (|4j) on a short time scale compared to the period of 
the laser field, we consider a static field approximation 
where we freeze the laser field at peak field amplitude: 



ns (x,p) 



£"02;. 



(5) 



Here we analyze this model for d = 1. It exhibits a 
Stark saddle [29] where the effective potential (Coulomb 
plus laser) is locally maximum and thus acts as a phys- 



ical barrier to ionization. In phase space, the saddle 
has zero momentum and its position IS a solution of 

2x*/ {xl + a^)^/^ ^Eo = 0. We denote = lis (^*, 0) 
its energy. Since the time dependence has been removed, 
the system has one degree of freedom and thus it is inte- 
grable. In phase space, the curve Hs = is a separatrix 
(homoclinic connection) between the nonionizing trajec- 
tories and the ionizing ones (see the outermost curve in 
Fig.Ej). 

Taking into account the time-dependent potential, the 
dynamics of Hamiltonian (j4]) is slightly complicated by 
the possibility of chaotic trajectories near the saddle 
point. In fact, from KAM type arguments, it is possi- 
ble to show the preservation of invariant tori (from the 
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FIG. 6. (color online) Phase space of Hamiltonian (O for d — 
1 and for a laser intensity / = 3 x 10^^ W cm~^ (surface). The 
separatrix associated with the Stark saddle 29] is represented 
as the outermost curve. For comparison, we display Poincare 
sections (stroboscopic plot at the maximum of the laser field, 
dots) for Hamiltonian (|4]) with d — 1. The inset shows a 
zoom of this Poincare section close to the Stark saddle. We 
represent the outermost chaotic layer [outermost dots (blue 
online)], the boundary torus [outermost light gray (magenta) 
curve in the bounded region], the elliptic island around the 
1:12 resonant periodic orbit [light gray (cyan) curves] and the 
hyperbolic 1:13 resonant periodic orbit. 



integrable dynamics described by Hs). We give an illus- 
tration of those tori in Fig. [6l Looking closer at the dy- 
namics, we notice some discrepancies between the static 
approximation and the time dependent dynamics. This 
is particularly the case close to the separatrix in Fig. [6l 
The bound region of Hamiltonian (j4]) corresponds to a 
part of phase space, close to the nucleus, filled up with 
KAM tori. The outermost tori corresponds to a bound- 
ary torus [outermost light gray (magenta online) curve 
in the inset of Fig. [6] outside which trajectories might 
be ionizing. This region is referred to as the unbound 
region, where the dynamics leads to ionization. Not all 
of the trajectories, though, ionize, as it can be seen for 
instance in a small elliptic island in the inset of Fig. [6l 
We notice two main dynamical features of the unbound 
region: Close to the boundary torus, we observe a thin 
chaotic layer [outermost dots (blue online) in the inset 
of Fig. [6]. Further away from the nucleus, the dynam- 
ics leads to rapid ionization (outermost white area in the 
figure). 

In the bound region of the dynamics defined by Hamil- 
tonian dll with one dimension, not all tori survive from 
the integrable case without the laser field. The weakest 
ones are the resonant tori. They are characterized by 
mvf = ncj, where is the free field rotational frequency 
of the electron around the nucleus and n,m G N*. They 
are later referred to as m:n resonances. According to the 
Poincare-Birkhoff theorem [30j , when a resonant torus is 



FIG. 7. (color online) Rotational frequency z/f of the elec- 
tron around the nucleus (dots) as a function of the maximum 
position Xm- For comparison, we display the theoretical pre- 
diction (dashed curve) obtained as described in the text. The 
comparison is done without laser excitation and a = 1. 



broken at least two periodic orbits are created instead. 
From numerical simulations, we find two or four resonant 
periodic orbits depending on the intensity. 

For the softening parameters a = 6 = 1, based on the 
observation that the rotational frequency depends almost 
linearly with energy, an analytical approximation for the 
action has been obtained in Ref. [22j]. The same formula 
can be generalized for other softening parameters, so that 
the approximate free field frequency z/f reads 



where 7, Sq and ujq are given by 



(6) 
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The maximum position Xm experienced by the electron is 
determined from E = + a^. As a consequence 

we obtain the expression of the rotational frequency z^f 
as a function of Xm- In Fig. [7] we compare this prediction 
with the actual frequency, computed numerically. Note 
the very good agreement between the curves, even rela- 
tively far away from the nucleus. 

Equation (|6]) combined with the resonance condition 
allows us to narrow down the regions of phase space in 
which to look for resonant periodic orbits. For the sim- 
plicity of the discussion we set the laser phase 0o = 7r/2, 
as most of our figures are drawn at peak field amplitude, 
such that t = corresponds to the peak field. In this 
case, for all l:n resonances (the electron oscillates ex- 
actly n times around the nucleus within one laser cycle), 
two resonant periodic orbits with zero initial momentum 
are present, one on each side of the nucleus. To predict 
their location, we use an adiabatic approximation: the 
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FIG. 8. (color online) Energy of the l:n resonant periodic 
orbits at peak field amplitude as a function of their position 
xq. Horizontal lines label the predicted position for the res- 
onant periodic orbits. Their actual position is depicted with 



dots. The parameters are /o 
a = 1. 



10^'' W cm" 



0.0584 and 



total energy for the inner electron Hin, in the bound re- 
gion, is almost constant in time. As a consequence, we 
replace the energy £ with Hin in Eq. (j6j) and the reso- 
nance condition becomes 



- EqXq — So 



(7) 



where xq is the position on the periodic orbit with zero 
momentum. Note that the implicit equation xq {Eq^uo^ n) 
has in general two solutions, one positive and the other 
negative. They correspond to two periodic orbits with 
zero initial momentum located on each side of the nu- 
cleus. In Fig.[8]we compare the results given by this pre- 
diction with the actual location of the resonant orbits. 
The very good agreement with the prediction is worth 
mentioning. For a given softening parameter a, there is 
a finite number of primary resonances (and thus a finite 
number of resonant l:n orbits): From Eq. ([7]) we deduce 
that the number of resonances nres is equal to 



2 1 



(8) 



where [-J is the floor function. For example, for a = 
1, we find nres = 24 primary resonances. This number 
decreases as a increases, i.e., as larger or equivalently 
more loosely bound atoms are modeled. 

For a given resonance index n, there are two periodic 
orbits with zero initial momentum, one on each side of 
the nucleus. We denote 0~ the one located on the same 
side as the Stark saddle (see the inset of Fig. [6]) and (9+ 
the other one. For the intensity range we have consid- 
ered, is found to be always elliptic while the stability 



of 0~ depends on the order of the resonance. For odd 
n, 0~ is hyperbolic for the whole range of intensities 
we have investigated. For even n, 0~ is elliptic for in- 
tensities where it is in the bound region and, at higher 
intensities, after being released in the unbound region, it 
bifurcates into a hyperbolic periodic orbit. This bifur- 
cation comes with the creation of an additional elliptic 
periodic orbit, through period doubling. For those reso- 
nances, in addition to the periodic orbits we find two 
additional hyperbolic resonant periodic orbits, which we 
denote O^. They are connected through a symmetry of 
the system, such that we identify the two and restrict the 
analysis to only one of them (see Fig. [51 the diamond in 
the upper right panel). 

As the laser intensity is increased, invariant tori are 
broken and the boundary torus between the bound and 
unbound region shrinks. This can also be seen in the 
static field approximation where this decrease is of the 
order of (which can be seen from the position of 

the saddle as given at the beginning of this section). At 
some point, the resonant periodic orbits are released to 
the unbound region. The dynamical influence of reso- 
nant periodic orbits depends strongly on their location 
and linear and nonlinear stability properties. Looking at 
Fig. [6] (inset), we note the elliptic island surrounding the 
1:12 resonant elliptic periodic orbit while the orbit is in 
the unbound region. Compared to Fig. O we note that 
this elliptic region corresponds to a part of phase space 
where no delayed ionization is detected. In the same fig- 
ure, the resonant orbit is hyperbolic. To analyze its 
nonlinear properties, we compute its stable and unsta- 
ble manifolds and display the resulting pictures in Fig. [5] 
(right panels). The manifolds are the pathways by which 
the electrons approach or leave the chaotic region. Note 
the strong similarities of the unstable manifold of 
with the Poincare sections of RESI trajectories in Fig. |H 
This similarity confirms the important role played by this 
unstable manifold (associated with a one-electron model) 
in the delayed double ionization process. 

Looking at Fig. O we see that the stable and unsta- 
ble manifolds associated with develop in the vicinity 
of the bound region and branches extend deep into the 
unbound region. Additionally, the width of the chaotic 
region as given by those manifolds explains why the tran- 
sition between the bound and unbound regions is so 
sharp [21]. A comparison between the left and right pan- 
els of Fig. [5] shows that the unstable manifold of reg- 
ulates the dynamics of the RESI trajectories. Even some 
of the very fine details of the unstable manifold are re- 
produced by this set of RESI trajectories. We note that 
some parts of the unstable manifold are missing: They 
correspond to the intersections between the stable and 
the unstable manifolds. The overlap between the stable 
and unstable manifolds of forms a "sticky" region [3l| 
that traps trajectories for some time before ionizing. 

For deeper insight into the shape of the unstable man- 
ifold of (9^, we represent a projection of this manifold in 
the plane (x, (j)) together with a projection of (black 
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FIG. 9. (color online) Projection of the unstable manifold 
of the periodic orbit O12 of Hamiltonian (|4| for / = 3 x 
cm ^ in the (x, 0) plane. Light gray corresponds 
to the manifold integrated for a long time while the dark 
gray focuses on the manifold integrated for a short time. The 
black curve is a projection of the periodic orbit O12 in the 
plane (x,0), and the black squares indicate the position of 
the Stark saddle [2^ at the maximum of the field. 



curve) in Fig. [9l It shows a central region near x = 
and two main branches which depart from this central 
region. It is worth noticing that these branches are lo- 
cated near the extrema of the electric field, i.e., when 
(p = 7t/2 and (j) = 37r/2. Since a delayed ionization tra- 
jectory follows this unstable manifold asymptotically, we 
see that this trajectory leaves the nucleus at times that 
correspond to the extrema of the electric field. This is 
consistent with what is already observed for RESI [12]. 
The unstable manifold has some other branches through 
which an electron can leave the nucleus and that do not 
correspond to extrema of the electric field. However, the 
probability of those ionization channels is much smaller. 

The discussion of the laser intensity considered in 
Figs. [5] and [6] can be generalized to arbitrary intensity. In 
this case, the hyperbolic resonant periodic orbit (s) in the 
unbound region {0~ or depending on the resonance 
index n) dynamically organize the chaotic layer surround- 
ing the bound region, through their stable and unstable 
manifolds. The unbound region hosts many resonant or- 
bits, associated with the accessible m:n resonances. How- 
ever, because of their short period (the same as the laser) 
and the short duration of the trapping for RESI (at most 
a couple of laser cycles) the primary resonances l:n are 
the ones which organize the RESI dynamics. The other 
periodic orbits in the vicinity of the bound region influ- 
ence only the fine details of the chaotic structure. As a 
signature of the influence of the l:n resonances, we rep- 
resent the normalized RESI statistics with respect to the 
double ionization yield as a function of the laser inten- 
sity for two-electron Hamiltonian ([1]) with d = 1 (lower 



FIG. 10. (color onhne) Normalized RESI yields for 460 nm 
wavelength laser (markers, left hand axis). Square (diamond) 
markers label one (two) dimensional one-electron simulations 
(normalized to the number of nonionized trajectories). Dot 
markers label two-dimensional two-electron simulations (nor- 
malized to the number of double ionizations) . The continuous 
and dashed curves represent the residue of the resonant peri- 
odic orbits On and of Hamiltonian Q (right hand axis) 
as a function of the laser intensity /. Continuous (dashed) 
curves refer to odd (even) l:n resonances. As intensity in- 
creases, the order of the resonance goes from n = 3 to n = 9 
as indicated on the curves. 

curve with square markers in Fig. [TQ|) . We clearly see 
some oscillations as a function of the intensity. In order 
to correlate these oscillations with some features of the 
l:n resonances, we represent the value of the Greene's 
residue [32] of 0~ and (9^ as a function of the intensity. 
These curves illustrate the strong correlation between the 
oscillations in the relative RESI yields with the successive 
releases of the resonant orbits in the unbound region. 

D. Dynamics of higher dimensional models 

In the adiabatic approximation, the one-electron 
model ([5]) is integrable for d = 1. A priori, this is not 
the case for two and three dimensions. In Fig. [11] we 
display the projection of the Poincare section for Hamil- 
tonian ([5]) with (i = 2, corresponding to the adiabatic 
approximation of Hamiltonian We see that the dy- 
namics is a mixture of regular and chaotic trajectories, 
with a large chaotic sea [light gray (cyan online) dots on 
the Poincare section] for energy values close to the Stark 
saddle. Since the laser is time independent, the separa- 
trix (red curve) acts as a barrier to ionization, and tra- 
jectories are confined inside the bound region (on each 
horizontal plane). Compared to Hamiltonian (j4]) with 
d = 1, the source of chaotic trajectories is very different: 
For Hamiltonian (j5j) with d = 2, it is the result of a non- 
linear coupling between the two spatial dimensions, while 
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FIG. 11. (color online) Poincare sections y = of Hamilto- 
nian ([5]) with d — 2 for a laser intensity / = 3 x 10^^ W cm~^ 
and a = 1 and energy ranking from the Stark saddle value 
(top) to 1-Ls — —2. The outermost (red online) curve corre- 
sponds to the one dimensional model with the same parame- 
ters (see Fig. [6]). 



for Hamiltonian (j4]) with (i = 1, it is the result of the in- 
teraction between the particle and the time-dependent 
field. It is expected that for the phase space of Hamil- 
tonian (|4]) with d = 2 these two sources of chaos both 
guide ionization. As a consequence, the dynamics allows 
for the possibility of an electron relatively deep inside the 
core region to diffuse chaotically and gain energy to be 
finally released after some delay. This phenomenon cor- 
responds to the mechanism we have referred to as chaotic 
diffusion when no collision or soft interaction is involved 
in the subsequent ionization. Compared to Fig. (TTJ the 
size of the chaotic region is amplified by the fact that 
additional tori are broken under the time dependence of 
the field. As a consequence, it becomes harder to define 
a bound and unbound region for the inner electron (see 
Fig. [12]), because of chaotic diffusion. 

As we did for the model @ with d = 1, the laser exci- 
tation can be considered as a perturbation of the free field 
dynamics for d = 2. Then, again, KAM type arguments 
show the persistence of invariant tori in the vicinity of 
the core. However, contrary to the one dimensional case, 
these tori are dimensionally too small to create barriers of 
transport in phase space and a chaotic diffusion can hap- 
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FIG. 12. (color online) Ionization time for Hamiltonian (|4|) 
with two dimensions for / = 3 x 10^^ W cm~^, 00 = 7r/2 and 
780 nm. Initial conditions are chosen at the energy of the 
Stark saddle 29] with px < and Py = 0. After a recollision 
(full line arrows), the inner electron can be directly ionized 
or promoted to an excited state that is ionized with a delay 
(dotted arrow). 



pen across them. To aid this view, and in a first (rough) 
approximation, one can imagine all the light gray (cyan 
online) dot regions in Fig. [11] to be connected, such that 
by moving from one layer to the next it is possible to 
cruise from one chaotic region to another even though 
there is a torus in between the two. More realistically, 
the manifolds of the resonant orbits, through the trans- 
verse direction(s) to the the field, develop relatively deep 
inside the core region. This allows for chaotic diffusion, 
which is drastically reduced for one dimensional models 
due to the existence of the boundary torus. Close to the 
nucleus however, because KAM tori form an intricate 
network, and by continuity of the dynamics, the diffu- 
sion is extremely slow (as stated in Nekhoroshev's theo- 
rem [33I). As a result, and since we consider very short 
laser pulses, in practice it is possible to identify a (finite 
time) bound and unbound regions with a chaotic layer 
surrounding the bound region. Compared to the one di- 
mensional case, this chaotic layer is more complicated, 
with branches developed deeper into the core region. 

In Fig. [To] we notice that the RESI yields given by 
Hamiltonian ([1]) are higher for d = 2 than for d = 1. 
We attribute these higher yields to the aforementioned 
chaotic dynamics coming from an increase of instability 
associated with additional resonances in higher dimen- 
sional models. It underscores once more the pivotal role 
played by the inner electron dynamics [Hamiltonian ^] 
in these delayed ionizations. 

For two dimensional models a direct consequence of 
the unstable direction for the resonant periodic orbits in 
the direction transverse to the laser field is the smooth- 
ing of the oscillations in Fig. [10] for 460 nm wavelength. 
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FIG. 13. (color online) Same figure as Fig. \T0\ for 780 nm 
wavelength laser. The upper panel corresponds to the stan- 
dard soft Coulomb potential while the lower panel corre- 
sponds to the modified potential [see Hamiltonian (|9))]. 



Increasing the wavelength to 780 nm gives birth to more 
resonances with a more dense tangle of unstable mani- 
folds so that the oscillations are completely washed out 
in two-dimensional models whether they are one- or two- 
electron models, as observed from our classical calcula- 
tions in Fig. [13] (upper panel). 

A second consequence of the unstable direction in the 
transverse plane is the increase of delayed double ioniza- 
tion yields as can be seen in Fig. [131 In the upper panel, 
the statistics for two dimensions (diamonds) is roughly 
one order of magnitude larger than it is for one dimension 
(squares). The same trend is confirmed with two electron 
models (compare Fig. [131 dots in the upper panel, with 
Fig. [151 filled in dots in the middle panel) with roughly 
one order of magnitude of difference between the one and 
two dimensional models. Looking at typical RESI trajec- 
tories we notice statistically longer delays for two dimen- 
sional models than for one. Again, we attribute these 
longer delays to the enhanced chaotic dynamics with two 
dimensions. 

Going to higher dimensions the difference between two 
and three dimensional models is not as drastic as between 
one and two dimensions because the two dimensions per- 
pendicular to the polarization direction are dynamically 
playing the same role due to the rotational symmetry of 
the system around the polarization direction. 

To illustrate the importance on the dynamics played 
by the direction transverse to the field, we consider a 
two-dimensional model where the Goulomb potential is 
modified in the ^/-direction. We consider the following 
modification of Hamiltonian ^ 



H(xi,X2,Pl,P2,t) = 



iPlI 



|P2| 



11^^X2 I 



Xi - X2 



62 



+ (xi + X2) • ej^Eof (t) sincjt. 



(9) 



where we Q is the matrix 




Tuning the parameter a enables one to modify the sta- 
bility of the relevant resonant periodic orbit in the trans- 
verse direction without changing its location and stability 
in the polarization direction. For numerical simulations, 
we set a = 2 and display the result in Fig. [131 (lower 
panel). For such a parameter, the resonant periodic or- 
bits become stable in the transverse direction. As a con- 
sequence of the stabilization of the transverse dynamics, 
the oscillations in the relative RESI yields are recovered 
(both with one and two electrons). Note also that the 
relative RESI yields are significantly reduced compared 
to Hamiltonian ([1]) which has a more chaotic dynamics. 



E. Significance of the target atom 

When changing the target atom, there has been shown 
an apparent strong sensitivity of RESI yields to the con- 
sidered atom [\ Isl, |7| . From the generic Hamiltonian ^ 
and its associated reduced models, these different atoms 
are modeled by varying the ground state energy and 
adjusting the softening parameters. Then, the analy- 
sis carried out in the previous sections can be applied 
to the specific atom. In particular, we focus on helium 
(a = 0.8), neon (a = 1) and argon (a = 1.5) and use 
the same electron-electron softening parameter for the 
three of them {b = 1). For instance, using Eq. (|8]) we 
find Tires = 33 primary resonances for He, rires = 24 for 
Ne and n^es = 13 for Ar so that the number of primary 
resonances l:n decreases with a as a~^/^. For all atoms, 
we find a similar organization of the delayed ionization 
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FIG. 14. (color online) Poincare section (dark dots, strobo- 
scopic plots at the maximum of the laser field, i.e., = 7r/2) of 
RESI trajectories of Hamiltonian (|4| with d = 1 for He^, Ne^ 
and Ar^ (from up to down) detected from a large set of ran- 
dom initial conditions for a laser intensity / — 10^^ W cm~^. 
On each panel, the color scale corresponds to laminar plots 
(distance, in logarithmic scale, of the electron from the core 
after 4.5 laser cycles [2l[) for the same Hamiltonian and the 
considered atom. 



dynamics by the resonant orbits (see Fig. [T4|). 

Previous studies have shown differences in the momen- 
tum spectra for various atoms 05 B • analysis sug- 
gests that this difference does not come from the phase 
space structures which are very similar for all atoms. To 
investigate the specific impact of the resonant orbits on 
delayed double ionization, we compare the RESI yields 
for each of the three atoms. The resulting picture is dis- 
played in Fig. [15] for Hamiltonian ^ and d = 1. As 
observed previously, in the nonsequential double ioniza- 
tion regime, the proportion of RESI exhibits oscillations. 
These oscillations are a clear signature of the resonant 
periodic orbits of the reduced Hamiltonian @ which or- 
ganize the chaotic layer around the bound region. Note 
the very good correlation between these oscillations and 
the residue curve of the l:n resonant orbits for the three 
atoms. To investigate the importance of the aligned ver- 
sus non-aligned model in delayed double ionization, in 
the middle panel we compare the RESI yields for Hamil- 
tonians (j2j) and (j3j). As expected, since they share the 
same inner electron reduced model, we observe the same 
trends between the two models where only the amplitude 
for low intensities is different. We also note that while 
similar oscillations are observed for the three atoms, the 
relative RESI yield is larger for Ne than it is it for both 
He and Ar, in agreement with what has been observed 
experimentally 
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FIG. 15. (color online) Same figure as Fig. \T0\ for 780 nm 
wavelength laser for Hamiltonian ([2|). In the panels we con- 
sider different atoms: He (upper, a — 0.8), Ne (middle, a — 1) 
and Ar (lower, a — 1.5). In the middle panel we also compare 
the relative yield obtained with the non-aligned model (|3)) 
(open circles). 



IV. CONCLUSION 

We have shown that delayed double ionization results 
from Hamiltonian chaos. This chaotic behavior results 
from well-identified resonances between the free field dy- 
namics and the laser excitation. Depending on the spe- 
cific question at hand, we have shown the limits of the 
reduced models: 

• When the preionized electron is far away from the 
nucleus, the dynamics of the remaining electron is 
accurately described by the inner electron reduced 
model (from which the resonances are identified). 
As a consequence, different two electron reduced 
models that share the same inner electron model 
[see Hamiltonians (j2j) and (j3j)] lead to qualitatively 
similar results for delayed double ionization. 

• While one-dimensional reduced models have been 
successful in describing direct impact double ioniza- 
tion, they give a partial description for the mech- 
anisms behind delayed double ionization. In this 
case, they exhibit the skeleton that regulates de- 
layed ionization through the l:n resonances. How- 
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ever, they miss the influence of the transverse di- 
mensions that enhance chaos. As a consequence, 
one dimensional models tend to significantly un- 
derestimate both the amount and the delays for 
delayed double ionization. 

• The mechanism of resonance is very robust 
against changes of parameters (laser frequency and 
species). In particular, for different atoms, the 
number of resonances varies, leading to a slight 
variation in normalized RESI yields. Nevertheless, 
in all cases, the resonances are correlated in the 
same way to the observed oscillations in the rela- 
tive RESI yields. 

All of these limits have been illustrated, e.g., on the plots 
of the normalized RESI yields versus intensity. In par- 



ticular, in some range of parameters, it has been shown 
that these curves exhibit oscillations which act as a clear 
signature of Hamiltonian chaos in these systems. 
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